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Abstract 

We propose a scheme to mimic the expanding cosmos in 1 + 2 dimensions 
in laboratory. Furthermore, we develop a general procedure to use nonlinear 
metamaterials to mimic D-brane and noncommutativity in string theory. 
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CsJ ■ 1. Introduction 

String theory is a developing field of modern theoretical physics which at- 
iy-j tempts to unify quantum mechanics and general relativity. The basic objects 

in string theory are one dimensional strings and high dimensional D-branes. 
To be self-consistent, the theory also requires supersymmetry, extra dimen- 
sions. The effective geometry on a D-brane becomes noncommutative when 
there is a nonzero background B field [l| . String theory has made tremendous 
theoretical progresses in the past 15 years, such as explaining microscopic ori- 
gin of the black hole entropy [2] in certain cases, the discovery of AdS/CFT 



[3| . However, string theory remains hard to be tested by experiments due to 
the extremely high energy scale. 

It will be of great interest if one can mimic some phenomena of string 
theory in laboratory. Fortunately, the dev elop ments of metamaterials and 
transformation optics 0, 0, @, 0, H S, E3, EH, S H Q E3 ma J be helpful for 



us to achieve this goal. Recently, metamaterials were used to design various 



interesting devices such as electromagnetic cloak |16|, [17[, perfect lens [18 



illusion devices [ 19l. |20| and so on. They can even be used to make an artificial 
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black hole |2ll l22l |23J and to mimic cosmos 24|, |25(. For example, in recent 



works [24|, |25], we find that the Casimir energy of the electromagnetic field 



in de Sitter space is proportional to the size of the horizon, the same form 



of the holographic dark energy j26| . We suggest to make metamaterials to 
mimic static de Sitter space in laboratory and measure the predicted Casimir 
energy 

In this paper we move on and propose to use metameterials to mimic 
some phenomena in cosmology and string theory. We design an approach to 
mimic a dynamic space-time in 1 + 2 dimensions by materials with constant 
permeability and varying permittivity, to avoid the complex issue of designing 



metameterials with synchronously varying permittivity and permeability 27 
Our proposal can be implemented by 3-level atoms or a nonlinear dielectric. 
Besides, we further suggest to use nonlinear metamaterials to mimic D-brane 
and noncommutative geometry in string theory. This topic is scientifically 
interesting, since it implies the possibility of testing and investigating some 
aspects of string theory at low energy scale in laboratory. Besides, the various 
properties of D-brane, noncommutativity and nonlinear electrodynamics will 
motivate us to design interesting optical devices. 

This paper is organized as follows. In Sec. II, we develop an executable 
experimental procedure to mimic expanding cosmos in 1 + 2 dimensions with 
abundant interesting effects such as redshift, the cosmic microwave back- 
ground (CMB) and the future event horizon. In Sec. Ill, we develop in 
theory a general procedure to use metamaterials to mimic gravity, D-brane, 
noncommutativity and so on. We conclude in Sec. IV. 

2. Mimic dynamic spacetime in 1+2 dimension 

In this section we discuss the issue of mimicking a 1 + 2 dimensional dy- 
namic space-time in metameterials. To do this let us consider the analogy 
between curved spacetime and dielectric medium for electromagnetic fields, 
as was explicitly studied in |10L Il6|. Especially, it is proven that light has the 
same behavior in dielectric medium and curved spacetime, under the con- 
ditions that the permittivity and the permeability are related to the metric 
by 



e a = £ _^LJ. g v i fi = -v-f-gv . (i) 

-goo soc^goo 

It should be stressed that this equation is valid only in 1+3 dimensions. 
This relation can be used to design metameterials and mimic specific kind 



of spacetime. For example, to mimic an expanding universe, the authors of 



271 ] proposed to vary the permeability by split-ring resonators, and change 
the permittivity by electrooptical modulation. 

Theoretically, one can make use of Eq.flT]) to mimic an arbitrary space- 
time. However, to mimic a dynamic spacetime, for instance, an expanding 
universe, one should simultaneously control the permeability and the permit- 
tivity to ensure them varying synchronously. Unfortunately, this is a hard 
task for experimental physics. To avoid this problem in this paper we just 
consider the case of mimicking a 1+2 dimensional spacetime. In this section 
we study two models, including a real 1+2 dimensional spacetime and a 1+2 
dimensional subspace of the 1+3 dimensional spacetime. In both cases we 
find that the permittivity and the permeability can be independent of each 
other, and the above problem is automatically resolved. Besides, it is natural 
and useful to consider the simplified case of realizing the 1+2 dimensional 
spacetime before turning to the much more complicated case of realizing 
a 1+3 dimensional spacetime. The metric of a homogeneous and isotropic 
universe in 1+2 dimensions takes the form 

ds 2 = g 00 c 2 dt 2 + g u dx 2 + g 2 2dy 2 , (2) 

with gu = g 2 2- Using the method developed in the next section we find that 
the dielectric with the following permeability and permittivity can be used 
to mimic this 1+2 dimensional metric 

/*=-^-, ^ = ^V(U = 1,2). (3) 

- £ o c goo —goo 

It is clear that the permittivity and the permeability are independent of 
each other. This remarkable property can be used to simplify experimental 
design. In fact, to mimic the expanding cosmos, one can keep the permeabil- 
ity as a constant and only vary the permittivity. We rewrite the metric in 
terms of the permittivity and the permeability from Eq.(j3]) 

ds 2 = ~^-c 2 dt 2 + ^JrW + dy 2 ) ■ (4) 

e [t) e[t) 

To realize the Friedmann- Robertson- Walker (FRW) metric 

ds 2 = -c 2 dr 2 + a 2 (T)(dx 2 + dy 2 ) , (5) 



we redefine the time coordinate as dr = e dt/e(t). Correspondingly we 
should express all the parameters in terms of this new time r rather than the 
laboratory time t. For example, fj,(r) = l/(eoc 2 ) and e(r) = eo/a 2 (r). 

Note that the Maxwell theory in 1+2 dimensions can be derived from that 
in 1+3 dimensions by imposing the self-consistent conditions E% = B 1 = 
B 2 = and d 3 Ei = d 3 E 2 = d 3 B 3 = 0. Hence, we can model this 1+2 
dimensional metric in dielectric with translation invariance along the z-axis 
and use transverse cross section perpendicular to the z-axis to mimic the 
two-dimensional space we are interested in. One may polarize photons to 
keep the magnetic field along z-axis, so that the magnetic field acts as a 
scalar on the cross section. Meanwhile, to confine them completely on the 
cross section, it is required that photons have no momentum along the axis. 
Therefore, one can make an analogy between the electromagnetic fields on 
the cross section and the counterparts in 1+2 dimensional spacetime 

Ei~Ei (i = 1,2), B 3 ~B . (6) 

Accordingly, one can also make an analogy between the relevant parameters 
of the dielectric and those in 1+2 dimensions 

e*i~e«(i,j = l,2), /i 33 ~/i. (7) 

Thus, to mimic the metric (JSJ), one only need to design metameterials with 
the property 

e» = e a- 2 (r) & (i,j = 1,2), jl 33 = l/(e c 2 ) . (8) 

We can also make the electric field parallel to the z-axis so that the 
electric field becomes a scalar on the cross section. In this case we can 
not mimic the Maxwell theory in 1+2 dimensions in which electric fields 
are vectors. Instead, we can use it to mimic the Maxwell theory in a 1+2 
dimensional subspace of the 1+3 dimensional spacetime. As in the above 
case, to be confined on the cross section, photons should have no momentum 
along the z-axis. When electric field is polarized along the axis, the nonzero 
components of the electromagnetic fields are E 3 and B^ (i,j = 1,2), and 
the relevant dielectric parameters are e 33 and ^ (i,j = 1,2). It should be 
stressed that if we consider only the 1+2 dimensional subspace, the irrelevant 
dielectric parameters such as e^ and u 33 need not take the exact form ([T|). In 
fact, they can take arbitrary forms and do not affect the polarized photons 



on the cross section. Setting g 33 = 1, one can derive the relevant dielectric 
parameters from Eq.flTJ) as 



£ 33 = £0V^ 5 ^ = _N/_0_« (ii = 12 ). (9) 

-goo so^goo 

We have the metric of the 1+2 dimensional subspace from Eq.Q 

d ^-W^) dt2 + Ww) idx2 + dy2) - (10) 

Next we set fi = l/(e c 2 ) and denote e 33 (t) by e(t), it follows that 

rfs 2 = _ c 2^2 + £W {dx 2 + rfj/2) (n) 

where £ is just the laboratory time. The scale factor, which signifies the 
scale of the universe, is a(t) = ^/e(£)/£ . Now we get the FRW metric in 
1+2 dimensional flat universe. One can mimic an expanding universe using 
metameterials with increasing e(t). 

In the rest of this section, we focus on the experimental design ffTTj) with 
the electric field polarized along the axis. Experimentally one has several 
alternative methods to vary the permittivity, e.g., by making use of the 
electrooptical effect and the model consisting of N identical 3-level atoms [28[ 
(the second method allows a comparably wider range of variation). Many 
interesting effects, such as the future event horizon, redshift and the varying 
CMB temperature, are expected in the experiment. 

If the permittivity increases fast enough, a future event horizon will arise 

rM ="Mfwr^f7m- (12) 

This is the largest distance that photons can travel. In the case of inflation, 
a(t) = e Ht (H is the Hubble constant), and rh(t) = cH^ 1 . Thus cH~ x is 
the limited distance light rays can propagate (the optical distance, not the 
coordinate distance). 

For light with frequency vq emitted at time to, its frequency at time t\ 
takes the form 29| 



i/i/i/o = a(t )/a(t x ) = Ve(W^i) • (13) 



If one increases e(t), the frequency decreases accordingly and light is red- 
shifted. The number density of photons is described by the Bose-Einsein 

distribution 

8tt u 2 du/c 3 ,„ ,. 

exp(hv/kj3T(t)) — 1 

Integrating this equation, one finds that the total photon number and the 
total energy iV oc T 3 and U oc T 4 , respectively. Thus the average energy of 
photons takes the form E = hv = U/N oc T. Notice that the temperature 



varies as T oc u oc l/a{t) oc 1/y e(£), which is exactly the same form of the 
CMB temperature. 

3. Mimic D-brane and Noncommutativity in String Theory 

In this section we demonstrate the feasibility of mimicking generalized 
electrodynamics by metamaterials. Let us start with the Maxwell equations 
in a dielectric medium 

d % E l = 0, e ijk djE k + d t B l = , (15) 

diD* = 0, e ij %H k - d t D { = . (16) 

where the electric displacement D l and magnetizing field Hi are related to 
the electric field Ei and the magnetic field B l through 

D t =e ij E jJ B i = /j ij H j . (17) 

For generalized electrodynamic theory in an arbitrary background spacetime, 
the Lagrangian takes the form LiF^^g^v). From the definition of the field 
strength F^ u = d^A u — d u A fM , one can derive the Bianchi identity 



d^F up + d„F m + dpF^ = . (18) 



By defining 



E % = F l0 , B* = ^e ijk F jk , (19) 

one can rewrite the Bianchi identity Eq. (ll8j) as Eq. (TT5|) . From the equations 
of motion (Eq.fjUJ) in the Appendix), one can derive Eq. (fT6]) as long as the 
electric displacement and magnetic field are defined by 

D = — M — ' Hi = d& — • (20) 



For details of the derivation please refer to the Appendix. Clearly, in the 
framework of the generalized electrodynamical theory, photons have the same 
behavior as those in the dielectric medium satisfying the constitutive rela- 
tions (120]) . due to the same expression of Maxwell equations in these two 
cases. This similarity implies a promising method to mimic generalized elec- 
trodynamics, D-brane and noncommutativity in string theory. 

Following our procedure, the well-known analogy ([Q) between curved 
spacetime and dielectric can be easily derived. For simplicity, here we take 
the case of 1+2 dimensions as an example. The Lagrangian of the electro- 
magnetic field in curved space is 

L(F) = -]e V=gF 2 = -]e Q ^g^ g ^F^F a p , (21) 



4 v a 4 

from which one can derive 

D i = ¥- = ¥~ = ^=V'^ , (22) 

dB dF 21 y/^ v ; 

Comparing the above equations with constitutive relations ( TT71) . one gets 

e« = ^El g v, » = ^-. (24) 

-#oo ~£oc z goo 

With this procedure at hand, we are in a position to investigate some 
interesting examples. In the following paragraphs, we study the mimicking 
of single Z) p -brane and noncommutativity in string theory. 

In string theory, the low-energy dynamics of single D p -brane is described 
by the Born-Infeld action 

S = -e T Dp J rf p+1 x^/-det(^ + 2vra'F^) , (25) 

where Tjj is the tension, p + 1 is the number of dimensions of the world- 
volume. From the above arguments, we can use an appropriate dielectric to 
mimic the low-energy dynamics of D p -bra.ne. Take p = 2 as an example, the 
corresponding dielectric has the following permittivity and permeability 

^oT^vrV 2 

fi fc . (26) 



4ttV 2 (c 2 £ 2 - E 2 ) 



As in the case of 1+2 dimensional expanding universe, one can mimic this D 2 - 
brane by dielectric with translation invariance along the z-axis and impose 
the self-consistent conditions E 3 = B 1 = B 2 = 0, d 3 Ei = d 3 E 2 = d 3 B 3 = 0. 

We encourage experimentalists to design such nonlinear metamaterials, 
they have great applications in mimicking string theory in laboratory. In 
such nonlinear metamaterials, a mass of novel phenomena are expected. For 
example, there is a upper limit of electric field intensity, the self-energy of the 
charge is finite, there is no birefringence and so on. To keep the permittivity 
and the permeability (|26|) real, it is clear that E < l/(2ira'). In general, when 
these parameters become complex, the system will exhibit an instability with 
abundant photons production and absoption. There is a correspondence of 
such instability in string theory. As we know, D-branes are objects where 
open strings can end, and charges carried by the two ends of a string are 
opposite. Thus, for an open string with both ends on the D-brane, an electric 
field will pull the two endpoints of string apart. If E is large enough, it will 
overcome the string tension to tear apart the string. It implies an instability 
similar to that in the dielectric. This instability of D-brane is hard to be 
observed in laboratory, but it is possible to mimic this effect by metamterials. 
Besides, metamaterials seem to have the unique capability to mimic such 
novel phenomenons. 

In string theory the effective geometry becomes noncommutative (NC) 
when there is a background 5-field [l| 

[/y^jT, (27) 

where x M are spacetime coordinates and 9^ v is a constant antisymmetric 
tensor. The generalized Yang-Mills action in NC space is 



S = --e l d i xF ta/ *F' M/ , (28) 

where the product * and field strength F^ are defined as 

f(x) * g(x) = exp(l^ A^_)/(a; + £) g ( x + C )| ?=c=0 , (29) 

F^ = d^K - d u A^ - i[A^ A„l . (30) 

In general, the NC Yang-Mill theory can be converted into a commutative 
but nonlinear counterpart by the Seiberg-Witten map. For example, when 

8 



Fn U is a constant, the Seiberg-Witten map is 



flV 



fll/ 



1 



he 



(31) 



where F a p is the ordinary field strength in commutative space. From the 
procedure developed in this section, one can mimic the NC Maxwell theory 
by metamaterials satisfying 



D' 



eoE* 



+ 



2Jfe (c 2 B 2 - E 2 )6 i0 



he ~ 



H 



1 + 2 V f£,0'° + eimn ^fe l ™B^Y 

e c 2 Bj 

1 + 1jfE0* + e lmn ^f6 lm B n ) 2 

(c 2 B 2 -E 2 )e m , 



e c- 



he 



1 + 2, 



'WEfl® + e lr 



ceo Qlm Qn\3 



(32) 



In NC geometry, Lorentz symmetry is broken. This can lead to birefrin- 
gence and direction dependent shift in the speed of light. For simplicity, we 
consider a constant electromagnetic field with small fluctuations, thus the 
Seiberg-Witten map ([3"Tj) for constant field strength is valid. Replace F^ u 
with fyw + F^, where f^ u denotes the constant background field and F^ u is 
the small fluctuation. The resulting Lagrangian reduces to 

£ „ nu .v £ 



4(1 + ^1^/^""" 



F, nJ F^ - 



- 2 Kx, u F^F^ 



where 



v af3"f5 



4(1 + v/p^ 



4 + x/fO^f, 



I fll/J 



it^f. 



(33) 



(34) 



he 



fW 



we keep only terms up to the second order in F^ u . The linear terms are 
absent due to the equations of motion of the background field. The equation 
of motion of the fluctuation is 



d v F^ + k^ pa d v F^ = . 



(35) 



For a plane wave F fiV (x) = F lll/ (p)e l P« a;a / R ) from the above equation one can 
get 

M^Ej = {-5 iJ p 2 - p>p k - 2k ifMUJ p f ,p u )E J = . (36) 

Following the general method, by requiring the determinant of M 4 - 7 vanishes, 
one can derive the dispersion relation of photons, then obtain the group 
velocity [30j 

v g = c\V0p°\ = c(l + p±<r). (37) 

to leading order in k^ upai where 

P = -\K, ° 2 = ^(M 2 - p\ (38) 

with 

As mentioned above, the velocity of light is direction dependent. Fur- 
thermore, since two solutions exist for E, it implies the birefringence effect, 
where these two components of light propagate independently with different 
velocities. 

To end this section, we list the main results in the leading order in the 
NC parameters and study an interesting optical phenomenon in such a NC 
space. It is likely that photons will spend less time traveling along polyline 
than along straight line. For simplicity, we focus on the case that all NC 
parameters vanish, other than 

# 03 = 9 > 0. To leading order of 6, the constitutive relations (1321) of 
metamaterials become, 

Hi = e c 2 (l + A x f^-6E z )B u 
V he 




D l = e (l + 4J^6E 3 )E l -2e (c 2 B 2 -E 2 )J^65 l 3 . (40) 



—* — * 

For simplicity, we consider only the case E — 0, B — (0, 0, B) for the back- 
ground field. Focus on the first order velocity (1371) . from Eq. (l34"l) and Eg. ( 1391) 
we learn that only p M in the leading order of 9 contribute, thus we can take 
]3 M = (1, sin<£>, 0, cosy?) in x-z plane as a good approximation {ip e [0, 7r/2] 
is the angle between the momentum of photons and background magnetic 
field). The velocity of photons then becomes, 




cos 2 <p) sin p . (41) 

10 



Figure 1: The polyline and straight line between O and A. 

We focus on the case v g = c + 2c^J~^0B\j / (1 + cos 2 if) sin 2 if below. It is 
clear that this velocity increases monotonically with if e [0,7r/2], thus the 
velocity assumes its minimum value v g = c when photons propagate along 
the direction of the magnetic field. It should be stressed that though the 
group velocity ( p4Tj) generally exceeds the speed of light c in vacuum, it does 
not mean that information can transfer faster than c. Various experiments 
have verified that it is possible for the group velocity of light in materials 
significantly exceed the speed of light in vacuum. However, in all these cases, 
the speed of information remains less than or equal to c [3j| . 

As depicted in Fig.l, the distance between these two points O and A is 
L. The time that photons will spend to cover this distance along the straight 
line is 

to = L/c . (42) 

Photons can also travel along another path, that is, the polyline depicted in 
Fig.l. In this case, photons are reflected to point A by a mirror, and the 
time it will take to travel along this polyline is 



h 



Lf cos ip 



c + 2c^f6B^(l + cos 2 if) sin 2 



H> 



(43) 



One can always find a nonempty interval for if to make sure t\ < to- For 
instance, set ^T^-QB = 0.1, one can find cti/L varies with if as depicted 
in Fig. 2. Hence, it takes less time for photons to move along the polyline 
than that along the straight line for if e (0,0.47). And the fastest path for 
photons to travel to point A is along the angle if = 0.24. Thus, the least 
time for photons to travel from point O to point A is t\ — 0.97L/c. 



11 




Figure 2: The time photons spend along the polyline with respect to the direction of 
propagation. 

It should be stressed that in such a NC space with a constant back- 
ground magnetic field, free photons still travel along straight line, but the 
fastest path for light to propagate between two points is generally a polyline. 
This interesting, counter-intuitive effect is due to Lorentz violation, which 
is believed to serve as signals of new physics coming from the Planck scale. 
Since the direct detection of Lorentz violation is difficult for the time being, 
it is of great significance if one can mimic and detect this effect in laboratory. 

4. Conclusion 

In this paper we discuss the issue of using metameterials to mimic some 
phenomena in cosmology and string theory. We propose a simple proce- 
dure to mimic a 1 + 2 dimensional expanding universe. We find that this 
can be implemented by materials with constant permeability and varying 
permittivity, so that the troublesome problem of designing metameterials 
with synchronously varying permittivity and permeability is automatically 
avoided. Moreover, we further suggest to use nonlinear metamaterials to 
mimic D-brane and noncommutative geometry in string theory. The permit- 
tivity and permeability of the corresponding metamaterials are calculated 
and provided in Sec. III. We see that it is possible to test and investigate 
some aspects of string theory at low energy scale in laboratory. However, it 
should be stressed that our proposals are as yet in theory, the metamateri- 
als which we suggest to mimic D-brane and noncommutative geometry are 
extremely difficult to make. We hope that with the progress in the field of 
metamaterials our proposals will be realized in the future. 
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Appendix 

In this appendix, we show that with definition ( 120]) . the Maxwell equa- 
tions in generalized electrodynamic theory and dielectric coincide. 

For a generalized electrodynamic theory, the Euler-Lagrange equation 

takes the form 

o dL dL 

d »w^7r3Ar°- (44) 

Note that 

ni_ dL dL dL dL 

one can easily obtain 

dT 

d - D ' - 9 -mA) = ° ' (46) 

Similarly, note that 

dL 1 dL ,_ 

Hk = ~dB^ = -r^m, ' (47) 

mnkTT _ _^,mi*, dL _ dL _ c dL 

k ~ 2 6 ^ k dF ZJ - C dF mn -2d(d n A m )> ^ 

from the Euler-Lagrange equation, one can derive 

^ H * - c 2 d -£r m - ~>^k = *>& = *°" • (49) 

Finally, with Eq. (j4U|) and Eg. ( 149]) . we get the same Maxwell equations f fT6l) as 
in dielectrics. It should be mentioned that a similar transformation-optical 



equivalence relations has been derived in 32 . 
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